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ABSTRACT 



With the recent results of helioseismology aboard 
SOHO, solar models are more and more constrained 
(Brun, Turck-Chieze and Morel 1998). New physi- 
cal processes, mainly connected to macroscopic mo- 
tions, must be introduced to understand these new 
observations. In this poster, we present solar mod- 
els including a turbulent pressure in the outer layers 
and mixing due to the tachoclinc (Spiegel and Zahn 
1992). 

Our results lead us to conclude that: 

- Turbulent pressure improves the absolute value of 
the acoustic mode frequencies (~ WfiHz at 4 mHz) 

- Mixing in a tachocline thickness of 0.05 ± O.O3i?0 
(Corbard et al. 1997) looks promising. 

Key words: sun: interior; sun: turbulent processes. 



1. THE SOLAR TACHOCLINE 



1.1. The Physical Description 

For most purposes, the Sun can be assumed in hy- 
drostatic and thermal equilibrium, neglecting the ef- 
fects of rotation and magnetic field. However models 
built with these simplifying assumptions do not agree 
with the helioseismic data, in particular with those 
obtained by the satellite SOHO, and it appears that 
macroscopic mixing processes must be taken into ac- 
count not only in the convection zone, but also in 
the radiative interior. Turbulent mixing may be in- 
troduced in the models by adding a turbulent term 
Dt in the diffusive part of the equation for the time 
evolution of the chemical abundance Xi. The equa- 
tion becomes: 



dXi d4Trpr 2 X l V i 



dt 



dm 



nuclear terms, 



(1) 



where the velocity Vi of species i with respect to the 
center of mass is: 



Vi = -^pr 2 {Di + D T 



d\nX t 
dm 



(2) 



The velocity Vi is the sum of one term which de- 
pends on the concentration gradient Di, and one 



which does not, Vi (Proffit and Michaud 1991). In 
this section, we address the mixing which may occur 
in the shear layer connecting the differential rota- 
tion in the convection zone with the solid rotation of 
the radiative interior. Spiegel and Zahn (1992) have 
given a physical interpretation of this tachocline by 
invoking anisotropic turbulence, with much stronger 
viscous transport in the horizontal than in the verti- 
cal direction. Such turbulence reduces the differential 
rotation and therefore inhibits the spread of the layer 
deep inside the radiative zone. The conservation of 
mass, momentum and entropy is given by: 



| + v.(pVO = o 



dV 



p I — + (V.V)V + 2flAV + flAr 



dt 



-VP-pV$ + V. || r 



V.V )S = V.(xVT) 



(3) 



(4) 



(5) 



where p is the density, P the pressure, T the tem- 
perature, V — (u, v, rQ sm9) is the local velocity in 

the rotating frame, (l the differential rotation, S the 
specific entropy, || r || the viscous stress tensor, $ 
gravitational potential. 

In order to solve the system of equations, some sim- 
plifying assumptions are made: 



• the flow field is axisymmetric V = V(r, 9, t) 

• anelastic approximation: dp/dt is negligible 

• advection is small compared to Coriolis acceler- 
ation 

• the tachocline is thin compared to the pressure 
scale height 

• viscous forces are small compared to Coriolis 
force. 



After separating each variable into its mean value on 
the sphere plus a perturbation, as T(r, t) + T(r, 9, t), 
the linearized form of the set of equations, in the 
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stationnary state (d/dt 
1992): 



0) is (Spiegel and Zahn 



P p T . 

— = — | equation ot state 

P p T 

IdP T . . 

— — — = g— hydrostatic equilibrium 
p or T 



-2£lrx£l 



1 dP 

pr dx 



geostrophic balance 



2S2x^- = p- 



vh{1-x 2 ) 



on 

dx 



dr dx 

diffusion and advection of angular momentum 

N 2 T d$> _ 1 

g pr 2 dx pc p r 2 



2 dT 



diffusion and advection of heat 



in spherical coordinates (r, x = cos 9), where ish is 
the horizontal turbulent viscosity (i>n >> vy) and 
defining a stream function W for the meridional flows 
by: 

9V . „ 

— , rpsmOv = — (6) 



r 2 pu 



We now project on horizontal cigenfunctions, the 

variables of the set of equations describing the 
tachocline: 

(P,f,u) = £(P„ 



= E* f Fn{x)da 



xtt = 



^n dK 



dx 



we find that the dominant term for the even eigen- 
function is n = 4. The final step is to connect the 
meridional velocity to the differential rotation. Ap- 
proximating the stream function by u 4 = ^4/ pr 2 and 
using the conservation of angular momentum, we get: 



_=0.5^p( M 4) ^ 



(7) 



with /i4 = 4.933. By introducing the variables ( = 
^(r bcz -r)/h and h = r^n/N) 1 / 2 ^/^) 1 / 4 the 
tachocline thickness, where rt, cz is the radius, N the 
Brunt- Vaisala frequency and K = x/p c p the radia- 
tive diffusivity at the base of the convective zone, we 
thus obtain this radial dependancy for 114: 



u A {r) = 0.5-f^/X4<5 4 exp(-C)cos(C). 

r bcz 



(8) 



The anisotropic diffusion invoked to stop the spread 
of the tachocline will also interfer with the advective 
transport of chemicals. Chaboyer and Zahn (1992) 
have shown that the result is a diffusive transport 
in the vertical direction, with an effective diffusivity 
given by: 



D T = 



D 



H 



E 



U 2 (r) 



n(n + l)(2n + 1) 



(9) 



where U n are the coefficients of the expansion of 
the vertical component of the velocity u in Legen- 
drc polynomials. The eigenfunctions F n above may 
be projected on these Legendre polynomials, and to 
a good approximation one has 



D T = 



r 



180 ' 



(10) 



Replacing U4 by (8), one reaches the following ex- 
pression for the vertical diffusivity 



Dt(C) = 



1 1 

1804 V 3 



h 



pt\ Olexp(-2C)cos 2 (C). 

(11) 



We shall treat the tachocline thickness h as an ad- 
justable parameter, roughly bounded between 0.03 
and 0.08 R Q (Corbard et al. 1997). With the new 
latitudinal dependence of £1 (Thompson et al. 1996), 
Vt bcz /2ir = 456 - 72x 2 - 42x 4 nHz, the coefficient 
Q 4 = 1.707 x 10~ 2 and fl/n = 0.9104. 



1.2. Results 



Starting from the reference model of Brim, Turck- 
Chieze and Morel (1998) built with the stellar evo- 
lution code CESAM (Morel 1997), which include the 
recent OPAL opacities and EoS, Alderberger et al. 
(1998) nuclear reaction rates and microscopic diffu- 
sion (see also Turck-Chieze et al. these proceeding), 
we introduce the coefficient Dt in our diffusion equa- 
tion with different values of h and N (fig. 1). The 
Brunt- Vaisala frequency N is taken as constant and 
h is chosen such as to yield a tachocline of w O.O5i? , 
as determined by (Corbard et al. 1997). In Table 1 
and figures 2,3,4,5 we present models including the 
coefficients shown in fig. 1. For one coefficient (Dti) 
we also test the influence on our results of the cali- 
bration of Z/X. 




Figure 1. Diffusive coefficients: Approximate value of 
the microscopic coefficient of A He (full line) , turbulent 
coefficient Dti (- - -), turbulent coefficient Dti (dot 
dash line) and turbulent coefficient Dt3 (dot dot dot dash 
line). 

Our standard model has a surface abundance of he- 
lium of 0.2425 in mass (see Table 1), which is a bit too 
low if we compare with the Basu and Antia (1995) 
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Table 1. Solar Models 

Parameters description: h: tachocline thickness, N: Brunt- Vaisala frequency, a: mixing length parameter, Yo, Zo, (Z/X)o: 
initial helium, heavy element and ratio heavy element on hydrogen, Y s , Zs, (Z/X)s: idem for surface compositions, Tb is 
the optical depth of the bottom of the atmosphere, Rf, zc , Tj, cz are the radius and temperature at the base of the convective 
zone, Y c , Z c , T c , p c : central helium, heavy element contents, central temperature and density. OPAL/A means that we 
use OPAL96 opacities and Alexander (1994) for low temperature. 



Parameters 


reference 


D T1 


D T1 


D T2 


Dt3 


Opacities 


OPAL/A 


OPAL/ A 


OPAL/A 


OPAL/A 


OPAL/ A 


Diffusion 


yes 


yes 


yes 


yes 


yes 


Age (Gyr) 


4.6 


4.6 


4.6 


4.6 


4.6 


h (r/Rn) 

\ 1 St). ! 




0.05 


0.05 


0.08 


0.05 


N (aHz) 




100 


100 


100 


10 


(Z/X) s 


fixed 


fixed 


free 


fixed 


fixed 


a 


1.768 


1.754 


1.763 


1.748 


1.759 


Y 


0.2721 


0.270 


0.2721 


0.2693 


0.2692 


Zo 


1.96 icr 2 


1.913 10~ 2 


1.96 10~ 2 


1.904 10" 2 


1.898 10~ 2 


(Z/X) 


0.0277 


0.0269 


0.0277 


0.0267 


0.0267 


Y 3 


0.2425 


0.2454 


0.2475 


0.2458 


0.247 


Z 3 


1.803 10~ 2 


1.796 10~ 2 


1.841 10~ 2 


1.79 10~ 2 


1.795 10~ 2 


(Z/X). 


0.0244 


0.0244 


0.0251 


0.0244 


0.0244 


Tb 


2 


2 


2 


2 


2 




0.713 


0.7145 


0.7134 


0.7145 


0.715 


T bzc x 10 6 (K) 


2.192 


2.181 


2.196 


2.181 


2.178 


Y c 


0.642 


0.639 


0.642 


0.638 


0.641 


z c 


2.095 10~ 2 


2.046 10~ 2 


2.096 10~ 2 


2.035 10~ 2 


2.03 10~ 2 


T c x 10 6 (K) 


15.71 


15.68 


15.71 


15.66 


15.67 


pc (g/cm 3 ) 


153.5 


153.3 


153.5 


153.0 


153.1 



value for the OPAL equation of state (Rogers, Swen- 
son and Iglesias 1996), Y s = 0.249 ± 0.003. When 
allowing for turbulent diffusion in the tachocline, the 
settling of helium is reduced by 16% for Dti up to 
25% with Dt3, leading to surface abundances Y s = 
0.2454 or Y s — 0.247 which are in better agreement 
with the helioseismic value. Note the smooth com- 
position profile in fig. 3, with its extended plateau 
below the convection zone. 
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Figure 2. Sound Square Speed Difference between 
GOLF+MDI data, the reference model (full line) and 
three models calibrated in Z/X — 0.0245 including a tur- 
bulent term: Dti ( ), Dt2 (dot dash line) and Dt3 

(dot dot dot dash line). 

We remark that the present Z/X value for a mixed 
model plays an important role in this study. If we 
do not calibrate the Z /X value of the mixed model, 



0.040 r 




o oon h ... i , ... i . ... i . ... i . ... ' 
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r/R s 

Figure 3. Hydrogen Composition for the reference model 
(full line) and three models with turbulent mixing due to 

the tachocline Dti ( ), Dti (dot dash line) and Dt3 

(dot dot dot dash line). 



starting with the heavy clement composition Zq of 
the reference model, we obtain different surface abun- 
dances (see for example column 2 and 3 of table 1 
for the coeficient Dti)- Also, when comparing the 
two sound speed profiles (figs. 4 and 5), we see that 
in both cases the bump around 0.7 Rq is pratically 
erased by the introduction of the tachocline mixing, 
but the calibrated model is modified along the whole 
structure and not only close to the bump. This could 
be understood by the fact that the tachocline mixing 
inhibits the elements settling and then less helium 
reaches the center, modifying the central conditions. 
The remaining part of the discrepancy between the 
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nuclear core and the convection zone may come from 
uncertainties on opacities coefficients or knwoledge 
of the element abundances (See Brun, Turck-Chieze 
and Morel f998) or from some other mixing. The 
low value of the tachocline thickness (h = 0.02) in- 
ferred by Elliot et al. (these proceedings) could be 
due to the turbulent coefficient used and also to the 
fact that they did not calibrate their model in Z/X. 
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Figure 4- Sound Square Speed Difference between 
GOLF+MDI data, the reference model (full line) and two 
models including the coefficient Dt\' calibrated in Z/X 
(- - -) and with a non calibrated Zq = Z^ d = 0.0f96 
(dot dot dot dash line). 




Figure 5. The effect of calibration on the sound square 
speed difference between the reference model and two mod- 
els including the coefficient Dti- calibrated in Z/X (full 
line) and with a non calibrated Zq = Z^ td — 0.0196 (- - 
-)■ 



2. EFFECT OF TURBULENT PRESSURE ON 
THE ACOUSTIC MODES 



We also introduce a turbulent pressure in the total 



pressure (i.e. P tot = P gas + P u 



turb 



where P gas is the 



sum of the gas+radiation+coulombian pressure). 



2.1. Implementation of the Turbulent Pressure 



The turbulent pressure is neglected in most solar 
models, in spite of the fact that it represents a 
substantial fraction of the total pressure in the up- 
permost part of the convection zone, which has an 



important weight in the determination of the fre- 
quency of the acoustic modes. By definition, the 
turbulent pressure is given by the horizontal aver- 
age Pturb —< >, v z being the vertical turbulent 
velocity, but here we link it to the convective veloc- 
ity v conv drawn from the Mixing Lenth Treatment 
(MLT): Pturb = Ppvl onv , where (3 takes into account 
the geometry of the turbulent motions. We follow 
Cox and Giuli (1968) who introduce a modified adi- 
abatic gradient . 



ad 



dlnP gaz 
dlnP 



(12) 



where V a£ ; is the adiabatic temperature gradient. 
Doing so we obtain for the convective efficiency T 
and for the velocity v conv a corrected expression: 



gSA 2 



V - Vl 



8A„ r + i 



(V-V ad ) (13) 



where, V is the temperature gradient of the medium, 
V' is the temperature gradient of the moving ele- 
ment, A is the mixing length, X p the pressure scale 
height, g the gravitational acceleration and S — 
— (dlnp/dlnT)p. After some algebraic manipula- 
tions, we deduce a new cubic equation for T: 



,r 3 + r 2 + r-A 2 [v ra(I -vJ = o 



(14) 



where V ra d is the gradient if all the energy is carried 
by radiation, A = CpKgS 1 / 2 p 5 / 2 A 2 / 12 vWZ^P 1 / 2 , 
0o = 9/4, c p the specific heat at constant pressure, 
a the radiation constant, c the light speed and k the 
opacity coefficient. The real gradient V is thus ob- 
tained by: 

, r(r + i) 

ad + (15) 



V 



We test different values for (3 from 0.5 to 1.5 (see 
fig. 7), and find that the best agreement with direct 
numerical simulations (see Rosenthal et al. (these 
proceedings)) are obtained for (3 = 1 (see fig. 6). For 
the other values, the turbulent contribution to the 
total pressure is either too low {Pturb ~ 8%) or too 

high {Pturb ~ 27%). 




Figure 6. Ratio of the total pressure Ptot = Pturb + Pgas 
over gas pressure for = 1. 
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2.2. Results 

By implementing the turbulent pressure in our solar 
model we obtain the following modifications: 

• a higher superadiabatic gradient at its peak near 
the surface 

• a reduction up to 10 \iHz of the 1=0 acoustic 
mode frequencies. 

This study confirms the important role played by 
the turbulent pressure in the upper layers of the so- 
lar convection zone (see Baturin and Mironova and 
Bohmer and Rudiger these proceedings). Rs intro- 
duction in the model removes a large part of the 
discrepancy between observed and predicted acous- 
tic frequencies, as shown in fig. 7. Again the best fit 
is obtained with (3 = 1. But this improvement can- 
not hide the fundamental shortcomings of the Mix- 
ing Length Treatment, in particular that of being a 
local theory which does not allow for convective pen- 
etration. Serious progress will only be achieved by 
realistic hydrodynamical modelization, as illustrated 
by the contributions of Rosenthal et al. (these pro- 
ceedings). 
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1.5 2.0 2.5 3.0 3.5 4.0 

Frequency (mHz) 

Figure 7. Frequency differences between GOLF 1=0 modes 
and theoretical model with turbulent pressure: (3 = 1.5 (- 
- -), (3 — 1 (dot) and without (full line). 
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